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ABSTRACT 

The long term (~ 1.5 years) X-ray enhancement and the accompanying in- 
frared enhancement hght curves of the anomalous X-ray pulsar IE 2259+58 fol- 
lowing the major bursting epoch can be accounted for by the relaxation of a fall 
back disk that has been pushed back by a gamma-ray flare. The required burst 
energy estimated from the results of our model fits is low enough for such a burst 
to have remained below the detection limits. We find that an irradiated disk 
model with a low irradiation efficiency is in good agreement with both X-ray 
and infrared data. Non-irradiated disk models also give a good fit to the X-ray 
light curve, but are not consistent with the infrared data for the first week of the 
enhancement. 

Subject headings: pulsars: individual (AXP IE 2259+58) — stars: neutron — 
X-rays: bursts — accretion, accretion disks 

1. Introduction 

Anomalous X-ray pulsars (AXPs) and soft gamma ray repeaters (SGRs) are neutron 
stars which are characterized by their persistent X-ray luminosities [L ~ 10'^^ — 10^^ erg s^^) 
well above their rotation powers (see Mereghetti 2002 and Hurley 2000 for reviews of AXPs 
and SGRs respectively). Another remarkable property of these sources is the clustering of 
their spin periods to a very narrow range (P ~ 5 - 12 s). There is no observational evidence 
for binary nature in these systems. They exhibit short 1 s), super-Eddington (<10^^ erg 
s^^) bursts. Burst repetition time scales vary from seconds to years. In about three decades 
of observations of these sources, three giant fiares were observed from the SGRs 0526-66 
(Mazets et al. 1979), 1900+14 (Hurley et al. 1999) and 1806-20 (Palmer et al. 2005). 
These giant fiares are characterized by an initial hard spike with a peak luminosity > 10"^^ 
erg s~^ which lasts a fraction of a second. The hard spike is followed by an oscillating tail 
that decays in a few minutes. Assuming isotropic emission, the fiuence of the entire giant 
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flare is >10^'' ergs (Hurley et al. 1999; Feroci et al. 2001; Mazets et al. 1999). Following 
the August 27*^ 1998 giant burst of SGR 1900+14, the persistent X-ray flux increased by a 
factor of ~ 700 to a peak luminosity ~ 10^^ erg s~^. The subsequent decay is a power law 
with index ~ 0.7 (Woods et al. 2001). The three other, less luminous X-ray enhancements 
of SGR 1900+14 were also preceded by gamma-ray flares. An X-ray enhancement was also 
observed from AXP IE 2259+58. Although no strong burst preceding the enhancement was 
reported, it is possible that such a burst could have been missed (Woods et al.2004). Out of 
six known AXPs including IE 2259+58, five sources were detected in the IR band (HuUeman 
et al. 2001; Israel et al. 2002; Wang & Chakrabarty 2002; Kaspi et al. 2003; Israel et al. 
2003; HuUeman et al. 2004; Israel et al. 2004). One of them, 4U 0142+61, was also observed 
in the optical (R) band (HuUeman et al. 2000, Dhillon et al. 2005). The optical flux from 
AXP 4U 0142+61 has a pulsed component with a pulsed fraction (~ 0.27) much higher than 
the observed pulsed fraction (~ 0.05) in X-rays (Kern & Martin 2002). Observations of 
transient AXP XTE J1810-197 separated by a few months showed a correlated 
decrease (by a factor ~ 2) in X-ray and IR luminosities (Rea, N et al. 2004, Israel 
et al 2004). IR radiation from this source was proposed to be radiated from a 
passive disk irradiated by the X-rays from the neutron steir Rea et al. (2004). 
In this model, the source of the X-rays is assumed to have magnetic origin, as 
proposed by the magnetar models, rather than accretion onto the neutron stcir. 
The X-ray enhancement of AXP IE 2259+58 was also accompanied by an enhancement in 
the IR luminosity (Kaspi et al. 2003). Long term (~ 1.5 years) IR and X-ray flux evolutions 
of this source exhibited a correlated decay behavior (Tam et al. 2004). 

The magnetar model (Duncan & Thompson 1992, Thompson & Duncan 1995) can 
explain the energies and the supcr-Eddington luminosities of the normal and the giant bursts 
of these sources. In these models, the source of the bursts is the huge magnetic energy release 
from inside the neutron stars. The persistent X-ray emission was suggested to originate from 
the magnetic field decay (Thompson & Duncan 1996). The spin down torque is provided by 
the magnetic dipole radiation. Magnetar models have no detailed predictions for the observed 
optical and infrared (IR) emission from AXPs and cannot explain the period clustering of 
AXPs and SGRs without assuming special conditions (Colpi, Geppert & Page 2000). 

Fall back disk models (Chatterjee, Hernquist, & Narayan 2000; Alpar 2001; Marsden et 
al. 2001) explain the period clustering of these sources using conventional dipole magnetic 
fields of young neutron stars with ~ 10^^ — 10^^ G on the surface of the neutron star within 
the time scales of the ages of their associated supernova remnants (Chatterjee & Hernquist 
2000; Ek§i & Alpar 2003). In fall back disk models, the source of the X-ray luminosity is 
the accretion onto the neutron star from the accretion disk. In earlier work (Ertan & Alpar 
2003), we showed, by means of numerical fits to data, that the X-ray enhancement of the 
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SGR 1900+14 following its giant flare can be reproduced by the relaxation of the inner disk 

subsequent to an initial push back caused by the flare. We estimated that part of the pushed 
back inner disk piles up and remains bound creating the post burst initial conditions of the 
disk, while the remaining part escapes from the system. The estimated energetics and the 
amount of escaping ejecta can account for the observed radio enhancement accompanying the 
X-ray enhancement of the SGR 1900+14 (Ertan & Cheng 2004b). An important limitation 
for the models is the observed optical and IR emission from the AXPs. It was shown that 
the observed pulsed component of AXP 4U 0142+61 can be explained by both the magnetar 
outer gap models and the disk-star dynamo models (Ertan & Cheng 2004a). Furthermore, 
unpulsed optical/IR emission from the AXP 4U 0142+61 is in agreement with the expecta- 
tions of a standard thin disk model (Ertan & Cheng 2004a). In the disk models, the inner 
radius of the disk is consistent with a dipole component of the magnetic field, with strength 
~ 10^^ G on the stellar surface. The origin of energetic bursts require magnetar strength 
surface fields. It could be the case that the burst involves surface magnetar fields in the 
higher multipole components, while a disk standing at inner radius set by a 10^^ G surface 
dipole field defines the rotation period of the star in a near-equilibrium state. It seems that 
such an hybrid model can explain almost all the observational behaviors of the AXPs and 
SGRs. 

The enhancement energetics of the SGR 1900+14 and the AXP IE 2259+58 
Eire different. While the SGR 1900+14 was observed only in X-rays, we now have 
contemporaneous X-ray and IR enhancement data for the AXP IE 2259+58. In 
this work, our aim is to answer: (i) Is it possible to explain the X-ray enhance- 
ment of the AXP IE 2259+58 by using a similar pushed back disk model applied 
to the SGR 1900+14? (ii) Is this model capable of reproducing both the X-ray 
and the IR data of the AXP IE 2259+58? In answering these questions, we 
test both irradiated and non-irradiated disk models. We estimate the required 
burst energy to build up post-burst initial inner disk conditions for the enhanced 
radiation and check whether such a burst could have remained below detection 
limits. 

The 2-10 keV X-ray fiux data of AXP IE 2259+586 were taken from Woods et al. 
(2004). They have monitored the source with the Rossi X-ray Timing Explorer / Propor- 
tional Counter Array, and determined the pulsed X-ray intensity of the source over a large 
time basehne. Using pointed XMM-Newton observations, they derived a conversion factor 
between the pulsed intensity and unabsorbed flux for this source. A detailed description of 
this conversion is given in Woods et al. (2004). 

The infrared measurements of IE 2259+586 were obtained from Kaspi et al. (2003), 
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Israel et al. (2003) and Tarn et al. (2004). The first two observations were made with the 
Near-Infrared Imager (NIRI) at the Gemini North telescope at about 3 and 10 days after 
the onset of the bursting activity. In these pointings, the Kg magnitudes of the source were 
found to be 20.36±0.15 and 21.14±0.21 respectively (Kaspi et al. 2003). The third infrared 
measurement, carried out with the Canada-France-Hawaii Telescope, gave the K' magnitude 
to be 21.31 ± 0.24 (Israel et al. 2003). The last two observations were also made with NIRI 
and resulted with magnitudes of 21.66±0.11 and 21.54±0.05 respectively. We converted 
these magnitudes into physical flux units using the method described in Tam et al. (2004). 
Details of the irradiated and non- irradiated disk models are given in §2. The results are 
discussed in §3. We summarize our conclusions in §4. 



2. The Numerical Model 

All four observed X-ray enhancements of the SGR 1900-1-14 were preceded by strong 
gamma-ray flares. No gamma-ray burst was reported before the onset of the X-ray enhance- 
ment of the AXP IE 2259-1-58. It is possible that such a burst could have remained below 
the detection limits. We assume that the X-ray enhancement of the AXP IE 2259-1-58 was 
preceded by a missed gamma-ray burst. We discuss the burst fiuence required to build up the 
post burst initial disk conditions and the observational limits in §3. In earlier work (Ertan & 
Alpar 2003), we showed that the X-ray enhancement evolution of the SGR 1900+14 following 
the August 27 giant flare can be accounted for by the relaxation of the inner disk matter that 
has been pushed back by the flare. We apply the same model to AXP IE 2259-1-58 assuming 
the physical mechanism giving rise to the X-ray enhancement is the same. Our aim is to 
explain the observed X-ray enhancement together with the accompanying IR enhancement 
data. In our model, the source of the X-rays is accretion onto the NS, while the IR emission 
originates from the accretion disk. We investigate the IR emission for both irradiated and 
non-irradiated disk models. 

In our numerical model, we take the model functional form of the post burst initial 
mass distribution to be the same as applied to SGR 1900-1-14 (Ertan &: Alpar 2003). The 
part of the inner disk matter pushed back by the flare energy which remains bound to the 
system is represented by a Gaussian T,{R,t = 0) = Umax exp — ( "^a^" ) ^ centered at the 
radius Rq. Added to this Gaussian is a surface density profllc in the form E = T,o{Ro/Ry 
which represents the extended outer disk. The inner disk radius R^^ where the inflowing 
disk matter is stopped by the magnetic pressure is kept constant. We numerically solve the 
disk diffusion equation for the surface density as described in detail in Ertan & Alpar (2003). 
We use the a prescription of the kinematic viscosity u — acsh (Shakura & Sunyaev 1973) 
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where cs = kTc/iimp is the local sound speed, Tc the local disk midplane temperature, 
h = Cs/f^j^ the pressure scale height of the disk and fl-^ the local Keplerian angular velocity 
of the disk. The viscosity parameter a = 0.08 and the outer disk radius -Rout ~ ^^^^ 
arc kept constant throughout the calculations. The chosen i?out large enough that it does 
not effect the results. 



2.1. X-ray Emission from the Neutron Stcir Surface 

For a thin disk, the total disk luminosity is L^^^-^ = GMM^^/2R^^, and most of this 
emission comes from the inner disk. Here, M^^ is the mass inflow rate arriving at the 
inner disk radius and M is the mass of the neutron star, which we take to be IAMq. 
The accretion luminosity from the NS surface, = GMM^/R^, determines the observed 
luminosity in the X-ray band. Assuming most of the X-ray flux from the source is emitted in 
the observation band (2-10 keV), we take the observed X-ray luminosity to represent the total 
luminosity of the neutron star. The evolution of the X-ray luminosity depends on both 
Mjjj(i) and its fraction accreted onto the neutron star / = M^/M-^^ where M* is the mass 
accretion rate onto the star. The X-ray luminosity can be written as = 2/(i?jj^/i?*)L(^jg]^. 
We compare the energy flux F^, ~ L^/^AikP) with the X-ray data for the model fits. We take 
the distance of the source d — b kpc (HuUeman et al. 2000). 



2.2. Infrared Emission from the Non-irradiated Disk 

For a thin, non-irradiated accretion disk, the source of the disk blackbody emission is 
the viscous dissipation inside the disk. The dissipation rate D per unit area sets the effective 
temperature and flux from both surfaces of the disk 

D = = 2aT^V (1) 

The surface density and the corresponding dissipation rate proflle at a given time is obtained 
by numerically solving the diffusion equation (see Ertan & Alpar 2003 for details). The 
effective temperatures are calculated for each radial grid point. We obtain the model IR 
luminosity by integrating the blackbody fluxes along the spatial grid points emitting in the 
observational IR band (K,, X = 2.15/im, width=0.31/im). The blackbody temperature 
which mainly emits in this band is kT^f^ ~ 0.2 eV. For the model fits, we relate 
the model IR luminosities to the observed fiux by Fm ~ Lircosi/ {And?) where i is the 
inclination angle between the normal of the disk plane and the line of sight of the observer. 
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The model puts an upper limit to the fraction / of the disk mass flow rate which is 
accreted onto the neutron star. This upper limit is obtained by setting cosi = 1. For the 
model given in Fig. 1, / = 0.03 and cosi = 0.9. The other model parameters are given in 
Table 1. It is seen in Fig. 1 that the non-irradiated disk model is consistent with the long 
term evolution of the IR data (~ 1.5 years) but not for the first two weeks. 



2.3. Infrcired Emission from the Irradiated Disk 

A likely reason for this discrepancy is the irradiation of the disk by the X-rays from the 
neutron star. For a disk irradiated by the X-rays from the neutron star, the irradiation flux 
through the disk surface must be added to the dissipative energy flux in Eq. 1. The X-ray 
irradiation flux can be written as 

^.^^ = <r = W 

where 

(Shakura & Sunyaev 1973), r] is the efficiency of the conversion of the rest mass energy 
into X-rays, e the X-ray albedo of the disk face, it! the radial coordinate and ifj^j. the 
pressure scale height of the disk which should be calculated including the effect of the X-ray 
irradiation (Dubus et al. 1999). The disk thickness to radial distance ratio H/R is roughly 
constant along the disk. The parameter C can vary in a large interval especially due to the 
uncertainty on the disk albedo rj. For a point irradiation source, estimates arc usually in the 
range 10"^ - 10"^ (Tuchman et al. 1990, de Jong at al. 1996, Dubus ct al. 1999). Using C 
as a free parameter, we calculate the effective temperature for the irradiated disk as 

<ff = f + ^irr- (4) 

Eq. 2 shows that Tj^.^. oc R-^^^, while T^g- oc i?^^/^ for a non-irradiated disk. For 
small radii viscous heating dominates over the X-ray heating. Beyond some 
critical radius Rc, depending on the irradiation strength, the main source of the 
emission from the disk surface is the reprocessed X-rays for an irradiated disk. 
The critical radius Rc can be estimated by equating the irradiation flux F^rr to 
the dissipative energy flux D/2 (Eq. 1). Assuming that the mass inflow rate at 
a given time is constant along the disk, we obtain 

3 CM* , /10-^\ „ 
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We follow the same method as that of the non-irradiated disk model to calculate the total 
IR luminosity in the observational IR (-ft's) band. In the vertical disk analyses, the X-rays are 
usually assumed to be absorbed in a thin layer of the disk surface. The irradiation modifies 
the vertical temperature profile and the effective temperature of the disk without significantly 
affecting the disk mid-plane temperatures that determine the kinematic viscosity (Dubus et 
al. 2001). Therefore, X-ray irradiation does not change the evolution of the disk mass 
transfer rates as long as the disk remains at the same viscosity state (see §3). The irradiated 
disk model cannot constrain the / parameter since there is the extra free parameter C . The 
model X-ray and IR flux curves given in Figures 2 & 3, were obtained with C — 2x 10~^ and 
/ = 0.26 for cosi — 0.76. When we take / near unity with cosi — 0.76, a similar model fit is 
obtained with C ~ 5 x 10~^. The annular disk section radiating in the Kg band is at a 
radial distance i? ~ 2 x 10^^ cm at the beginning of the enhancement, and shrinks 
slowly to i? ~ 5 X 10^° cm along the decay phase, remaining in the irradiation 
dominated part of the disk. Figures 2 & 3 show that this particular irradiated disk 
model is in agreement with both the X-ray and the IR data. 



3. Results and Discussion 



The model X-ray and IR light curves are presented for the non- irradiated (Fig. 1) and 
irradiated (Figures 2 & 3) disk models. The corresponding model parameters are given in 
Table 1. The differences between the initial surface density distributions of the best fitting 
models are their amplitudes and a small change in the power- law index p of the extended disk 
profiles. In the irradiated disk model, since the effective temperatures are mainly determined 
by the irradiation flux rather than the viscous flux, the amplitude of the initial surface density 
distribution is less than that of the non-irradiated disk model (by a factor ~ 5). The time 
evolution of the mass accretion rate M onto the neutron star which is observed as the 
X-ray light curve is similar for all different models of the disk. For a given initial mass 
distribution, the evolution of the mass flow rate M through the disk is the same for both 
models, as the irradiation has no significant effect on the disk midplane temperatures. In 
the hydrogen disk models employed to account for the observed outbursts of dwarf novae 
and soft X-ray transients, the disk has two states. The viscosity parameter a is large (small) 
when the effective temperature of the disk is high (low) compared to a critical temperature 
reflecting the ionization of hydrogen (~ 10^ K). A local disk instabihty can lead to an 
outburst (quiescence) if the resultant heating (cooling) wave can take all or most of the disk 
to the high (low) viscosity state. The inclusion of the X-ray irradiation in models modifies 
significantly the critical accretion rates at which these instabilities occur. Nevertheless, for a 
given stable state, the effect of irradiation on disk mass fiow evolution is small. The critical 
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mass flow rates and the critical temperatures for disk instabilities depend on the 
exact composition of the disk material which determines the opacities in the 
disk (Menou et al. 2002). They use carbon - oxygen composition to discuss 
disk instabilities. This composition is uncertain for AXP and SGR disks. In all 
our calculations, we use the electron scattering opacity as we previously applied 
for SGR 1900+14. During our simulations we kept the viscosity parameter constant at 
the value a. — 0.08 which is characteristic of the hot states in hydrogen disk models. We 
also tested disk models with both high and low viscosity states occurring (bimodal) going 
through the hydrogen ionization temperature. We used a parameters similar to those of 
hydrogen disk models (ct/jot ~ 0.1 and otcou ~ 0.01 — 0.05). For a critical effective 
temperature ~ 10^ K, both X-ray and IR model curves deviate from the observed trend of 
the light curves significantly following the instability. The X-ray enhancement light curve of 
the SGR 1900+14 can also be reproduced by using a constant a parameter a ~ 0.1 (Ertan & 
Alpar 2003). Secular decay of the light curves imply that in the accretion regime 
indicated by the observed X-ray luminosities of the AXPs and SGRs, fall-back 
disks do not experience a global disk instability. 

The rise and early decay phase (the first few days of enhancement) of the X-ray luminos- 
ity are determined by the release of the Gaussian surface density distribution representing 
the post burst pile up. In this phase, X-ray model curve evolution depends mainly on the 
amount of mass included in the Gaussian and the distance between the Gaussian center (i?o) 
and the disk's inner radius, rather than the detailed Gaussian parameters. The two model 
light curves seen in Fig. 2 cire for two different initial Gaussian surface density 
distributions with the same extended (outer) disk profile. The constant inner 
disk radius was taken near the Alfven radius for a dipole magnetic field with 
strength ~ 10^^ Gauss on the stellar surface and Mjj^ ~ 10^^ g s^. Similar fits 
can also be obtained for larger R^^i ^o- The model fits do not constrain 

the dipole magnetic field strength. For the first few weeks, the decay phase does not 
strongly depend on the chosen power-law index p of the extended model disk surface den- 
sity profile, provided it remains roughly around its best fit value for the overall evolution 
curve. About a few weeks after the onset of the enhancement, light curve evolution becomes 
sensitive to the power index p. For a steady disk, p is expected to be around 3/4, taking 
the run of the midplane temperatures as i?"^/^. This value of p is close to our best fit p 
parameters for both irradiated and non- irradiated disk models (0.65 and 0.70). In Fig. 
3, we illustrate three irradiated disk model light curves for different p values 
(0.60,0.65 and 0.70), keeping the Gaussian parameters constant. Irradiated and 
non-irradiated models eire successful in explaining the long term X-ray enhance- 
ment data. The non-irradiated disk model can reproduce the IR data except 
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for the first two weeks with a low upper limit / < 0.03 for the accretion rate to 
disk mass flow rate ratio /. The irradiated disk model gives a better fit to the 
overall IR data without constraining /, due to the uncertainty in the efficiency of 
the X-ray irradiation for the AXP and SGR disks which is likely to be different 
from the LMXB disks. For LMXBs optical to X-ray luminosity ratios are much 
higher than expected optical luminosity from the intrinsic dissipation of the disk 
and the X-ray luminosity from the neutron star surface. While this indicates 
that the optical fiux of LMXBs cire from an X-ray irradiated disk, theoretical 
calculations of the disk thickness profile show that irradiation must be indirect 
because of the self-screening of the disk (See e.g. Dubus et al. 1999). Most 
likely source of the indirect emission is a hot scattering corona above and below 
the central disk region. The hot corona can be fed by the thermally unstable, 
optically thin surface layers of the hot innermost disk. Is it plausible to expect 
that the SGR and AXP disks have the same X-ray irradiation efficiency as that 
of LMXBs? This is not likely because the stronger magnetic fields of SGRs and 
AXPs cut the disk at radii two to three orders of magnitude IcU-ger than the 
inner disk radii of LMXBs. For the irradiated disk model, we initially fix the 
peu-ameter / = M * /^in obtain a good fit to the X-ray data points 

first. Then, we look for the best fit to the IR light curve tracing the irradiation 
strength C. The irradiated disk model curves presented in Figures 2 and 3 are 
obtained with / ~ 1/4, cosi ~ 3/4 and C = 2 x 10""^. We see that fits with similcO" 
quality can be obtained for fCcosi ~ 4 x 10~^. 

The gamma-ray burst energy required to build up the initial conditions of our models 
can be estimated by assuming the post burst pile up represented by Gaussian surface density 
distribution was accumulated near the inner disk radius just before the burst. The mass 6M 
of the Gaussian distributions is ~ 8 x 10^° g and ~ 1.6 x 10^*^ g for the non-irradiated and 
irradiated models respectively. The amount of energy needed to push these masses from 
the inner disk radius to Ro obtained from the best model fits are 2.6 x 10^^ ergs and 
5.2 X 10^^ ergs. For an isotropic burst, the fraction of the burst energy absorbed by the 
disk is 27ri?(2/i)/(47ri?^) — h/R which is nearly constant along the disk and about a few 
xlO~^ for the accretion regime of the AXP IE 2259+58. Then, the required total burst 
energies are about 10^° ergs for the non-irradiated disk model and 2 x 10^^ ergs for the 
irradiated disk model. If the burst took place before the RXTE observations presented here, 
a burst duration more than about 20 seconds for the non-irradiated disk is enough for the 
burst luminosity to remain under the detection limit of the Ulysses gamma-ray detectors 
which were active before the RXTE observations (Woods et al. 2004). For the burst energy 
expected from the irradiated disk model, the burst duration must be longer than only ~ 5 
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seconds for the burst not to be detectable by the Ulysses gamma-ray detectors. 

4. Conclusion 

We have shown by means of numerical fits to data that the long term (~ 1.5 years) 
contemporaneous X-ray and IR energy flux data of the AXP IE 2259-1-58 can be accounted 
for by the evolution of the disk after inner disk matter is initially pushed back by a burst. 
The estimated burst energy from the initial conditions of both irradiated and non-irradiated 
disk models arc small enough to remain below the detection limits of the Ulysses gamma-ray 
detectors operating before the enhanced X-ray observations. 

Both irradiated and non-irradiated disk models can reproduce the X-ray enhancement 
data. The irradiated disk model is also in agreement with the overall IR hght curve ac- 
companying the X-ray enhancement. The non- irradiated disk model gives a slower decay 
than that indicated by the first two IR data points corresponding to about three days and 
ten days after the onset of the X-ray enhancement. Nevertheless, it is consistent with the 
subsequent 1.5 year long term IR data. Since the early decay phase could depend on the 
details of the initial disk conditions we can not exclude the non-irradiated disk model. 

It is not possible to compare the active and passive disk models on the basis 
of the IR (K band) data (for a given X-ray light curve). Since the emission 
in this band comes from the irradiation dominated outer part of the disk, the 
same IR radiation would be expected to be observed in either model. However, 
expectations of the two models are different for the shorter wavelength optical 
emission (R, V, B bands) and their relation to the X-ray flux. For an active disk 
with a mass inflow rate ~ 10^^ — 10^^ g s~^, significant fraction of the emission in 
these optical bands comes from the intrinsic dissipation near the innermost disk 
radius. Dissipation gives rise to certain amount of mass inflow, and thereby to 
accretion since this mass inflow is near the co-rotation radius for relevant mass 
inflow rates. Our model thus expects a correlation also between the dissipation 
dominated optical (R,V,B bands ) luminosity and the X-ray luminosity. (Note 
that an inner disk emitting in the optical indicates that the dipole magnetic field 
is ~ 10^^ G for relevant disk mass inflow rates). Furthermore, any cheiracteristic 
VEiriation in the dissipation dominated part of the optical luminosity is expected 
to be followed in X-rays with a delay chciracterized by the viscous time scale 
of the inner disk (hours days in persistent states, much shorter in the early 
enhancement phase due to large pressure gradients). For AXPs, in the persistent 
phase, our model can estimate the dissipation dominated optical flux using X-ray 
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, IR and inter-stelleir extinction information, which requires a detailed source by 
source examination of the AXPs (work in progess). For AXP IE 2259+58, we 
estimate that the luminosities in and R bands are nearly the same, and their 
ratios to the X-ray luminosity is about 10 in its persistent state. At present, we 
cannot use this point to distinguish between our model and passive disk models 
like that of Rea et al.(2004), because only upper limits exist for the R and V 
band luminosity of AXP 2259+58. However, this remains as a distinguishing 
prediction of our model. 
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Fig. 1. — X-ray (plus signs) and IR (triangles) energy flux data of the enhancement phase 
of IE 2259+58. The solid lines are the non-irradiated disk model curves. 
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Time (days since 52443.65 MJD) 

Fig. 2. — X-ray (plus signs) and IR (triangles) energy flux data of the enhancement phase of 
IE 2259+58, and irradiated disk model curves. Dashed and solid curves are for two different 
initial Gaussian surface density distributions (see Table 1 for the model parameters) . 
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Time (days since 52443.65 MJD) 

Fig. 3. — X-ray (plus signs) and IR (triangles) energy flux data of the enhancement phase 
of IE 2259+58. The irradiated model curves correspond to different power law indices p of 
the initial extended outer disk surface density distribution. The other disk parameters are 
the same for all three curves. 
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Table 1: Model parameters for the flux evolution presented in Figures 1 & 2. The three model 
curves in Fig. 3 have the same parameters as given in the first column of the irradiated 
disk model, but for different p values (0.60, 0.65, 0.70). See §3 for the explanation of the 
parameters here. 





Non-irradiated disk Model 


Irradiated Disk Model 


T^max (g cm 2) 


7.5 X 10^ 


1.5 X 10^ 


3.5 X 10^ 


Gaussian width (cm) 


2.0 X 10' 


2.0 X 10' 


1.1 X 10' 


So/Smax 


0.6 


0.6 


0.3 


i?o (cm) 


1.0 X lO'^ 


1.0 X lO'^ 


1.0 X lO'' 


^in (cm) 


8.3 X 10^ 


8.5 X 10^ 


8.5 X 10^ 


P 


0.7 


0.65 


0.65 


f 


0.03 


0.27 


0.27 


cosi 


0.9 


0.76 


0.76 



